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In the last years different studies have revealed the usefulness of a microcanonical analysis of 
finite systems when dealing with phase transitions. In this approach the quantities of interest 
are exclusively expressed as derivatives of the entropy S = In D. where f2 is the density of states. 
Obviously, the density of states has to be known with very high accuracy for this kind of analysis. 
Important progress has been achieved recently in the computation of the density of states of classical 
systems, as new types of algorithms have been developed. Here we extend one of these methods, 
originally formulated for systems with discrete degrees of freedom, to systems with continuous 
degrees of freedoms. As an application we compute the density of states of the three-dimensional 
XY model and demonstrate that critical quantities can directly be determined from the density of 
states of finite systems in cases where the degrees of freedom take continuous values. 
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I. INTRODUCTION 

In the microcanonical treatment of a finite system the 
main quantity of interest is the specific entropy s(e,m) 
as a function of the energy e and of the order parameter 
m Q. The latter is given by the magnetization in the 
case of a magnet. Recent investigations of finite classi- 
cal systems with discrete degrees of freedom undergoing a 
phase transition in the infinite volume limit have revealed 
that the microcanonically defined spontaneous magneti- 
zation m sp (e) is zero for energies larger than a certain 
energy e c and rises steeply with a power law behaviour 
for energies smaller than e c 0,0] Interestingly, the cor- 
responding susceptibility, being directly related to the 
curvature of the entropy surface, diverges at e = e c . It 
has to be noted that the exponents governing the power 
law behaviour of the different quantities in the vicinity 
of e c take on the classical mean field values for all sys- 
tem sizes smaller than infinite 0]. The true nonclas- 
sical exponents of the infinite system can, however, be 
determined from a microcanonical scaling analysis @. 
Whether such a behaviour with a sharp onset of the or- 
der parameter and a diverging susceptibility is termed 
a "phase transition in the microcanonical ensemble" or 
not, is a semantic question and, as such, of lesser im- 
portance than the question if there really exists a point 
with a true divergence already for a finite system. This 
has frequently been challenged with the reasoning that 
for a discrete system as e.g. the Ising model the argu- 
ments of the entropy assume only discrete values. For 
an Ising system with nearest neighbour interactions on 
a d-dimensional hypercubic lattice with linear extension 
L these values are ek — Ek/N and mi = Mi/N with 
Ek = 4:kJ and Mi = 21, where k and I are integers, J is 
the coupling constant and N = L d is the number of spins 
in the system. Only in the limit of infinite system size do 
the ratios of differences become equal to the derivatives 



(with respect to e and m) needed for the calculation of 
the susceptibility and other physical quantities. 

In order to overcome these serious objections against 
the microcanonical way of analyzing critical phenomena, 
we have decided to determine the entropy of a classical 
spin system where the spins take on continuous values. 
Of course, there the entropy s(e, m) is a continuous func- 
tion of its arguments. The system chosen is the XY 
model in d = 3 dimensions which undergoes a second 
order phase transition. 

In recent years different numerical methods have been 
proposed for the computation of the density of states of 
classical models @,l*u1u1l|]- Here we consider a highly 
efficient algorithm which has been applied to several dis- 
crete models as e.g. the two- and three-dimensional Ising 
models 0], the three-states Potts model 0,0], the vector 
Potts model with four states [1] or the voter model [Tlj . 
For the computation of the density of states of systems 
with continuous degrees of freedoms we have to modify 
this algorithm. 

The outline of the paper is the following. In Section 
2 we generalize the numerical method presented in 0] to 
models with continuous variables. This method yields 
the transition variables which permit the construction of 
the entropy surface. It is important to note that tran- 
sition variables can also be obtained when using other 
algorithms as for example the Wang-Landau method 
Therefore our method for deriving the entropy from these 
quantities can be applied very generally. Data obtained 
in this way for the three-dimensional XY model are an- 
alyzed microcanonically in Section 3. This enables us to 
determine critical exponents in systems with continuous 
variables directly from the density of states, either by ex- 
trapolating effective exponents or from a microcanonical 
finite-size scaling ansatz. Section 4 gives our conclusions. 



2 



lcro state 



C 



-si -< 



c 

_ _ -c 



(V) 



(K) 



FIG. 1: The density of states at different macrostates [i, v 
and K (schematic). For a reversible mechanism the number 
of connections C = C(/i — > v) = C(f — > fi) between two 
macrostates is the same for the forward and backward tran- 
sitions. 



II. THE COMPUTATION OF THE 
MICROCANONICAL ENTROPY 

A. Transition variables 

The algorithm developed in the following allows the de- 
termination of the density of states (DOS) (and there- 
fore also of the microcanonical entropy) of classical spin 
systems. In an extension of the method introduced in 
Q, the computation of the microcanonical entropy is 
performed in the case where the DOS is a function of 
continuous variables. The method is exemplified for the 
three-dimensional XY model with the classical Hamilto- 
nian 



u = - ^2 St- Sj 



(1) 



where the spin Si , characterizing the lattice point i of a 
simple cubic lattice, is a two-dimensional vector lying on 
the unit circle. The sum in Q is over nearest neighbour 
bonds. 

The consideration underlying our method is relatively 
simple. Assume at first a classical discrete spin sys- 
tem. A macrostate of the system may for convenience 
be denoted by \i. One might think of the energy and 
magnetization mi that characterize a macrostate or level 
jU = (efe,m;). In general, a huge number of microstates 
belongs to a macrostate /x. This number is the degener- 
acy Z^ of the level [i. In the course of the simulation a 
reversible mechanism takes the system from a microstate 
in a level \i to a different one which belongs to another 
level v. This step is repeated many times. Of course, 
the same mechanism has to be applied at every update. 
Starting from one microstate, a number N of new mi- 
crostates can be generated. The number N depends on 
the mechanism and on the model under consideration. 

When the mechanism operates on all microstates of 
the level [i then MZ^ new microstates can be gener- 



ated. A number C(/x — > v) of these belong to the level v. 
The quantity of specific interest is w(fJ, — ► v) = C(fi — > 
v)/AfZ[j, which is the probability of arriving at anyone of 
the microstates belonging to the level v when the start- 
ing point was one of the microstates of level ^. As the 
mechanism is reversible the numbers of forward and back- 
ward transitions between the levels \x and v are the same: 
C(fi — > v) = C{y — > /x), see Figure This then yields 
the expression 



w(y — > fi) 



(2) 



which allows the determination of the degeneracies 
from the transition probabilities w. Eq. (j2j, which is 
sometimes called Broad Histogram Equation @|, is al- 
ready complete for discrete spin models like the Ising or 
Potts models. 

Let us now focus on spin models with continuous spin 
variables. As an example we discuss the XY model [12j . 
but the generalization to other models (as for example the 
Heisenberg model) is straightforward. In our case both 
the energy and the magnetization have continuous values 
and the magnetization is a two-dimensional vector. As 
the energy function (QJ is invariant under global spin ro- 
tations, it is sufficient for many investigations to consider 
only the modulus m of the magnetization. This also con- 
siderably reduces the amount of resources needed for the 
numerical simulation of the model. For the purpose of 
performing this simulation we discretize the energy and 
the magnetization, the discretized values being denoted 
by e and to, and call Se respectively Sm the width of the 
discretization. We thereby allocate all microstates with 
energies e between e — Se/2 and e + Se/2 and modulus m 
of the magnetizations between m — 8m/ 2 and m + 8m/ 2 
to the same macrostate denoted by fx = (e, to). As the 
magnetisation m = (toi,TO2) is a two-dimensional vector 
the "volume" 



V(m) = 2nm8e8m 



(3) 



in (e,m)-space is not a constant. The DOS tt(e,m) is 
derived from 



fl(e,m) = Z(e,m) / V(m) 



(4) 



where Z(e, rh) = Z^. This expression also holds for other 
models beyond the XY model, where one only has to 
replace the volume © by the appropriate expression. 
It is important to note that during the simulation the 
spin variables as well as the energy and the magneti- 
zation of the system adopt continuous values. The dis- 
cretisation only concerns the quantities depending on the 
macrostates, e. g. the DOS and the transition probabili- 
ties. 

As already mentioned we perform the simulation with 
a reversible mechanism. For the XY model we use single 
spin rotations of randomly selected spins and a random 
rotation angle — tt < ip < ir. The following description of 
the algorithm is very general and also applies to discrete 
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spin models. Suppose that the system is in a macrostate 
denoted by e and rh or equivalently by [i. In the next 
step it is attempted to bring the system into another 
state v by the use of the mechanism described above. 
We increase the number of attempts B{\i) to leave the 
macrostate ^ by one count. At the same time we add 
one count to T(fi — > v) which is the number of attempted 
transitions from [i to v. In a long run the ratio t(/i — * 
v) := T(n — > v)j B{n) finally approximates the transition 
probability u>(/i — ► v). It is the quantity — ► v) which 
we call transition variable. During the simulation B(fi) 
and T(fi — > v) are updated at every attempted step. The 
probability of acceptance of the transition from jj, to v is 
chosen to be 



V 



t(/J, -> v) ' 



(5) 



This choice of p leads to an equal number of attempts 
B([i) to leave any macrostate /i, if sufficiently many up- 
dates of the system are performed. This is due to the 
fact, that now the probabilities for an actually executed 
transition for both the forward and backward directions 
are equal. It is worth mentioning that the probability 
of acceptance changes during the simulation and that it 
approaches its asymptotic value 



p = mm 



w(v — > p) 



w{n — > v) 



(C) 



for very long runs. 



B. Construction of the entropy surface 

Suppose that a macrostate fi = (e, m) whose degener- 
acy has not yet been determined can be reached from 
several macrostates v and suppose that the values of the 
Z v are known, then with the knowledge of the transition 
probabilities one may calculate lnZ^ starting from one 
of the macrostates with v = v': 

w(v' — * fi) 

In Zn = In — + In Z,,> 

— > v') 

and one would obtain the same result lnZ M for each of 
the states v as a starting point. However, the simula- 
tion only yields the transition variables t(fi — > v) which 
are estimates for the transition probabilities w{^jl — > v). 
Consequently, as these estimates are subject to stochas- 
tic fluctuations, we end up with different values for lnZ M 
when starting from different macrostates v. 

To take this into account we propose to estimate In Z^ 
by the weighted average 



In Z,, 



E 



Pi* 



In 



t{p -> v) 



InZ,. 



(7) 



where the sum is over all macrostates v from which the 
macrostate /i can be reached. The weights p^ are given 



by 



with 



Tjp -> v)Tjv -» \i) 
T(/i-»i/)+T(i/-»/i) 



(8) 



(9) 



which follows from the gaussian error propagation. An- 
other possible choice for the <7 M „ is given by 



g^v = min (T(/x -> v);T{v -> /i)) . 



(10) 



We checked that both choices JjjJ and <ff3~0f> lead to mi- 
crocanonical entropies which, within the errors, are the 
same. 

Up to now we have tacitly assumed that the degenera- 
cies Z v entering the sum in Eq. (JJ were known exactly. 
This is of course not the case, the Z v being also affected 
by statistical errors. To reduce the errors in the entropy 
values we propose the following iterative scheme. At the 
start we attribute an arbitrary value hiZ Va = Co to a 
randomly chosen macrostate vq. Note that the choice of 
Co does of course not affect our analysis as microcanoni- 
cally defined physical quantities only involve derivatives 
of the microcanonical entropy Q. Furthermore, if one 
would use the generated entropy in a canonical analy- 
sis, the constant Co would drop out when taking canon- 
ical averages of any observable. Having chosen an initial 
macrostate, we perform a random walk in the parameter 
space (in our case the energy-magnetization space) and 
attribute to any visited macrostate /i the degeneracy Z ^ 
obtained by evaluating Eq. 0. Of course, in doing so, 
only macrostates which have in the past been visited at 
least once can contribute to the sum in Eq. Q. This is 
repeated until every macrostate has often been visited, 
typically a few thousand times. When In Z^ has been ob- 
tained for all the macrostates fi = (e,m), Eq. (QJ yields 
the microcanonical entropy as (setting ks = 1) 



Sin) = \nZ lx -lnV(m). 



(11) 



It is of advantage not to integrate directly over the 
whole parameter space, but to restrict the random walk 
at the beginning to a small area around the starting point 
vq, as shown in Figure [2 The area at the start must be 
large enough to encompass all the macrostates reachable 
from the macrostate Vq. Typically after a few thousands 
visits of every macrostate inside the initial area, this inte- 
gration area is increased. This change of integration area 
is performed periodically. At some stage we no longer 
update the center of the integration area, see Figure EJ 

The presented method has the advantage that all tran- 
sition variables computed in the simulation are used for 
the construction of the entropy. Furthermore, the it- 
erative scheme just described ensures that the error in 
the final result for the entropy is minimized. It must be 
noted that this method is of course also applicable for 
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FIG. 2: The integration area at the beginning (left figure) and 
after several extensions (right figure). For this sketch we sup- 
posed that the initial point is only reachabel from microstates 
inside the grey area (left). At later stages, the middle of the 
integration area is not changed anymore (right). 



other models, including models with discrete spins, as 
only the correct volume V has to be inserted into Eq. 
(Till . Finally, we mention that for larger system sizes it 
is often of advantage to compute the entropy of stripes 
restricted in energy direction @, [3. The adaption of 
the presented method to this restricted geometry in the 
parameter space is straightforward. 



III. ANALYSIS OF THE ENTROPY 

The procedure outlined in the previous Section has 
been applied to systems of different linear extension L 
up to L = 25. For most of the runs the size of the chan- 
nels has been chosen to be unity on the extensive scale, 
i.e. Se = 6m = l/N, though different discretizations 
have also been tried with no noticeable effect on the re- 
sults. For L = 10 this is illustrated in Figure where 
we show the computed entropy as a function of m for 
two different fixed values of the energy and two different 
discretizations. The data obtained for the different dis- 
cretizations perfectly agree within errors. Note that this 
is an important point as the independence of the results 
on the bin size demonstrates that we are indeed accessing 
continuous properties. 

Whereas for small systems the entropy s(e,m) can be 
computed in a single run, for larger systems the increas- 
ing number of channels in energy and in magnetization 
direction makes the computation of the entropy very time 
consuming. One can then restrict the determination of 
the transition variables to narrow stripes restricted in 
the energy direction, which permits a distribution of the 
computation on different CPUs. The width of our stripes 
was usually 25 units on the extensive energy scale, which 
corresponds to twice the maximum energy increment in a 
single move. On the intensive scale the stripes are rather 
narrow (« 3 • 10" 3 ) for a system with L = 20. Therefore 
the resulting entropies can be averaged over the 25 en- 
ergy channels in order to obtain s(e,m) at the center e 
of the stripe. 

Having the entropy or, equivalently, the density of 
states at our disposal we can use these data in different 
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FIG. 3: Computed entropy for the three-dimensional XY 
model with N = 10 3 spins as function of the modulus of 
the magnetization for two different energies and two different 
discretizations: 8e = 8m = l/N (full symbols) and Se — 
8m = 2/N (open symbols). The data have been obtained 
by simulating two different stripes centered at the energies 
e = —0.65 and e = —1.5. Only selected points are shown by 
symbols. Note that for e = —0.65 the entropy is maximal 
for m = 0, whereas for e = —1.5 the maximum is located 
at m = m 3p ,L(e) > 0. For illustrative reasons the data have 
been shifted vertically by adding a constant. 



ways. Of course, we can compute the canonical partition 
function from the density of states and then obtain ther- 
mal averages of the different quantities of interest as the 
energy, the susceptibility and so on. Here, we take a dif- 
ferent route and directly investigate the microcanonical 
entropy s(e,m) itself. Recent investigations of discrete 
spin systems 0, 0- EiJ E^l have shown that critical expo- 
nents governing the power law behaviour of quantities in 
the vicinity of a critical point can be determined reliably 
by analyzing directly the density of states of finite sys- 
tems. Here we extend this analysis to phase transitions 
taking place in classical spin systems with continuous de- 
grees of freedom This kind of approach also poses a 
severe test to our method of computing the entropy, as a 
microcanonical analysis requires data of extremely high 
quality due to the absence of the smoothening effect of 
the Boltzmann weights. 

Coming back to our numerical data we observe that the 
entropy s(e, to), as a function of the modulus of the mag- 
netization, shows one maximum at to = for e > e C: L 
and one maximum at m = m sp x{e) > for e < e c see 
Figure^! For a given energy e the DOS as function of m 
exhibits an extremely sharp maximum such that the over- 
whelming majority of the accessible states belongs to the 
value of to where s has a maximum. Therefore m sp ^{e) 
is identified with the spontaneous magnetization of the 
finite system HUES- Figure shows the spontaneous 
magnetization for three system sizes. In order to pro- 
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FIG. 4: Spontaneous magnetization of the three-dimensional 
XY model as function of the energy per spin for three differ- 
ent sizes. Note the sharp onset of the spontaneous magneti- 
zation at the size-dependent energy e Cj _L. Error bars are much 
smaller than the sizes of the symbols. Lines are obtained by 
fitting the data close to e c ,£ to a square root 1131 . 



duce this plot the maximum has been determined for a 
series of stripes centered at different values of e. Notice 
the sharp onset of m sP) .£,(e) at e Cj L, a behaviour also ob- 
served for discrete spin models displaying a continuous 
phase transition in the infinite volume limit @, Q0| 
This onset is accompanied by a divergence of the radial 
part of the microcanonically defined susceptibility [T^ | 
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(12) 



This divergence is due to the fact that the curvature in 
the radial direction of m (i.e. the denominator in l|12(l ) 
vanishes at e Cj £. 

Exploiting the fact that close to e c ,£ the spontaneous 
magnetization rises with a square root singularity, e Cy L 
can be determined with high precision from a fit of the 
nonzero values of m sPj i(e) by the expression 



isp,h{e) = A L (e c , L - e) 



1/2 



(13) 



where Al and e Cj L are adjustable parameters. The re- 
sulting values of e Cj i, are plotted in Figure^las a function 
of 1/L together with the fit function 



e c ,L 



B(l/L)"' 



(14) 



with the adjustable parameters e c ,oo, B and v e . The fit 
yields the values e C;00 = —0.965(10) and v e = 0.666(5) 
which, considering the relatively small systems we have 
simulated, are in good agreement with the best values 
found in the literature, namely e c — —0.9884 and 
v = 0.67155(27) 0. 



PIG. 5: Pseudo-critical energy e c ,L, defined by the sharp 
onset of the spontaneous magnetization, as function of the 
inverse linear extension 1/L. The full line results from a fit of 
Eq. l[l4j to the data points, with e c ,oo = -0.965, B = 4.534 
and Ve = 0.666. Error bars are comparable to the symbol size. 



We briefly pause to recall that in the entropy formal- 
ism considered here critical exponents are not always 
identical to the better known thermal critical exponents. 
Indeed, in cases where the specific heat is diverging at 
the critical point (implying that the thermal critical ex- 
ponent a governing this divergence is positive, a > 0) 
critical exponents x e appearing in the microcanonical 
analysis are related to their canonical counterparts x by 
x e = x/{l-a) 0. This is not so in cases where the spe- 
cific heat displays a cusp singularity with a < 0. There 
x t = x. For the three-dimensional XY model we have 
a —0.015 fl9l | and consequently critical exponents ob- 
tained in the microcanonical analysis are identical to the 
thermal critical exponents. We can therefore savely drop 
the index e from now on. 

For all finite system sizes the expansion \1'M . which is 
valid in the vicinity of e c ,L, yields the classical critical 
exponent /3 C) £ = 1/2 Q|, but the range of validity of 
(FLUl shrinks to zero with L — > oo. Outside of this range 
TO sp.L(e) differs very little from m sp (e) = m sp _ 00 {e) for 
not too small system sizes. The critical expansion for the 
infinite system is given by 

Ae? 



m sp (e) 



(15) 



where e 



t^- and e g = —3 is the ground state energy 
per degree of freedom. Eq. iflfill is stricly speaking only 
valid in the limit e — ► 0. One may analyze the sponta- 
neous magnetization bylooking at the energy dependent 
logarithmic derivative j3| 



Peff{e) 



dim 



<sp 



dln(e) 



(16) 



which approaches the true critical exponent (3 in the limit 
limit e — > 0. Plots which show this approach for the two- 
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FIG. 6: Effective exponents derived from the spontaneous 
magnetization, see Eq. Ilfip . as function of e for different sys- 
tem sizes. The dashed line extrapolates the data for L = 20 
with 0.15 < e < 0.50 to zero. 



FIG. 8: Microcanonical finite-size scaling plot for the three- 
dimensional XY model. The values j3/v = 0.53(1) and 1/v = 
1.47(2) result from the best data collapse. Error bars are 
comparable to the sizes of the symbols. 




leads to the scaling ansatz 

I? 1 " m sp , L ((e c , L -e))~W (c (e c , L - e) L 1 '^ (17) 

where W is a universal scaling function characterizing the 
given universality class and C is a nonuniversal constant. 
It is worth noting that in this approach the knowledge of 
the value of the critical energy e c of the infinite system is 
not needed. As shown in @Eq. H2J) yields reliable values 
of the involved critical exponents for discrete models. We 
show in Fig.|Hlthat this is also the case for the continuous 
model discussed here. Plotting lP/ v m sPl L against (e c ,L — 
e)l}l v ', we obtain the values (3/v = 0.53(1) and 1/v = 
1.47(2) from the best data collapse, in good agreement 
with the expected values 0.5189 and 1.4891. 



FIG. 7: Variation of (3 £i l obtained from extrapolating linearly 
the effective exponents fi e ff as function of L. The dashed line 
indicates the literature value /3 = 0.3485 [Tfl| . 



and three-dimensional Ising systems can be found in [3j. 
Figure shows /3 e //(e) of the XY model under study for 
several lattice sizes. All the graphs, of course, precipitate 
to zero on approaching e = because m sp ^{e) is finite 
at e = e c , i.e. for e = 0, but the region of constant slope 
around e — 0.5 can be used to extrapolate the data to 
e = 0. Figured shows the result of this extrapolation for 
five values of L. These values quickly come close to the 
expected value (3 = 0.3485(2) 

Microcanonical finite-size scaling, as suggested in is 
another way of determining critical exponents. This ap- 
proach takes advantage of the existence of a well-defined 
transition point e c ^ in finite microcanonical systems and 



IV. CONCLUSIONS 

In recent years impressive progress has been achieved 
in the numerical computation of the DOS of finite classi- 
cal systems, as new and highly efficient simulation meth- 
ods have been proposed. One of the aims of this work is 
the generalization of the method presented in 0] to sys- 
tems with continuous degrees of freedom. This method 
relies on the powerful concept of transition variables @|. 
In our implementation the actual values of the transition 
variables are used for the acceptance rates of a proposed 
move, and this during the entire run. As a consequence 
all sites of the parameter space are visited with the same 
rate. 

With the final values of the transition variables we can 
construct the density of states, using the complete in- 
formation contained in these variables. We have formu- 
lated this approach in the present paper for magnetic 
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systems, where the natural variables are given by energy 
and magnetization, but a generalization to other situa- 
tions is straightforward. 

With the knowledge of the density of states we can 
analyze the system under investigation in different ways. 
One common approach is to compute the partition func- 
tion and then proceed with a canonical analysis. The 
good quality of our data, however, also enables us to per- 
form a microcanonical analysis where derivatives of the 
entropy with respect to energy and magnetization play 
a predominant role. The usefulness of the microcanon- 
ical analysis in the study of phase transitions has been 
revealed in many recent studies Q, 0, 0, 0, |2iJ, HU| • 

The theory of phase transitions is usually formulated 
within the frame of the canonical ensemble, which is 
based on the partition function or on its logarithm, the 
free energy. All the thermodynamic functions, as e.g. the 
susceptibility, are calculated from derivatives of the free 
energy with respect to the temperature and the applied 
fields. For all finite system sizes these derivatives remain 
finite, but in the thermodynamic limit they may diverge 
at the critical point, if a continuous phase transition oc- 
curs in the system. 

In the microcanonical ensemble the analysis of phase 
transitions is based on the DOS or its logarithm, the 
entropy Q, 0, EH ■ It has been demonstrated for several 
discrete spin models that the abrupt onset of the order 
parameter at e c and the divergence of the susceptibility 
occurs already for systems consisting of a rather modest 
number of spins, albeit with the classical values of the 
critical exponents 



We have demonstrated in this work that similar fea- 
tures are observed in systems with continuous spins: In 
a finite system of linear extension L the curvature of the 
entropy surface in the magnetization direction changes 
sign at the point e = e Cj j,. It is there where the order 
parameter sets in and where the susceptibility diverges 
in continuous as well as in discrete systems. Therefore 
for the XY model, where e and m take continuous val- 
ues, the entropy s(e,m) shows the same behaviour as for 
discrete spin models where it is defined only at discrete 
points in the parameter space. 

Of course, numerical studies must be performed at dis- 
crete values of the variables. The distinguishing fact is 
that for continuous models these values can be chosen 
freely, such that derivatives can, in principle, be calcu- 
lated from the ratios of arbitrarily small differences. The 
reality of a numerical study, however, sets the limits. 

As a final application we have demonstrated that in 
systems with continuous degrees of freedom critical ex- 
ponents can in principle be determined directly from the 
density of states, along the same lines as in discrete mod- 
els 

HUGHE!- In particular, we obtain very good esti- 
mates of v and f3 with regard to the modest system sizes 
used in this study. 
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